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Abstract: We use a field theoretic generalization of the Wigner-Weisskopf method to study 
the stability of the Bunch-Davies vacuum state for a massless, conformally coupled interact- 
ing test field in de Sitter space. We find that in A^ 4 theory the vacuum does not decay, while 
in non-conformally invariant models, the vacuum decays as a consequence of a vacuum wave 
function renormalization that depends singularly on (conformal) time and is proportional to 
the spatial volume. In a particular regularization scheme the vacuum wave function renor- 
malization is the same as in Minkowski spacetime, but in terms of the physical volume, which 
leads to an interpretation of the decay. A simple example of the impact of vacuum decay 
upon a non-gaussian correlation is discussed. Single particle excitations also decay into two 
particle states, leading to particle production that hastens the exiting of modes from the de 
Sitter horizon resulting in the production of entangled superhorizon pairs with a population 
consistent with unitary evolution. We find a non-perturbative, self-consistent "screening" 
mechanism that shuts off vacuum decay asymptotically, leading to a stationary vacuum state 
in a manner not unlike the approach to a fixed point in the space of states. 
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1 Introduction 

The stability of de Sitter space has always been an interesting problem to ponder, whether in 
the context of inflation or of the cosmological constant problem. Early studies [1-3] revealed 
that de Sitter space time contains infrared instabilities and profuse particle production in 
interacting field theories [4]. Particle production in a de Sitter background has been argued to 
provide a "screening" mechanism that leads to relaxation of the cosmological constant [5-7] 
much like the production of particle-antiparticle pairs in a constant electric field. Cosmological 
expansion modifies the energy-uncertainty relation allowing "virtual" excitations to persist 
longer, leading to remarkable phenomena, which is stronger in de Sitter space time [8]. 
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Polyakov[9] has recently argued that the Bunch-Davics vacuum for interacting test scalars 
is unstable to vacuum decay. This process induces runaway particle production, thus render- 
ing the in and out vacua inequivalent. Polyakov then argues that the effect of this process 
is to counteract the cosmological constant, not unlike what happens in the Schwinger mech- 
anism, leaving an non-inflationary FRW universe in its wake. While Polyakov's arguments 
are compelling, a firm confirmation of this effect with a sound and systematic calculational 
framework is still lacking. 

The arguments in ref. [9] are bolstered to a certain extent by a variety of calculations 
of various decay rates for particle states in de Sitter space. In refs. [10, 11] as well as in 
ref. [12], the rate for a field to decay into its own quanta was calculated, while the vacuum 
decay rate for a massless conformally coupled scalar with a quartic potential was computed 
in refs. [13-15]. These latter authors use an S-matrix approach to compute the decay rate 
for |0) — > \k±, . . . ,k^) and argue that the lack of energy conservation allows this decay to 
proceed. However, it seems to us that these calculations need to be taken with a grain of salt. 
At the outset, the use of S-matrix amplitudes presupposes the existence of asymptotic states 
connected by an S-matrix that may not even exist in de Sitter space[16]. More generally 
though, if the background is time dependent and does not turn off in the asymptotic regions, 
physical quantities will be time dependent and must be tracked as such. 

Turning to the explicit decay rate calculations, we know that for a massless conformally 
coupled with a A0 4 potential, we can perform the conformal rescaling to x = a (v) c l ) i where 
a(rj) is the scale factor in conformal time i] to convert this to the case of a massless interacting 
field in Minkowski space. But the Minkowski vacuum is perturbatively stable. How can this 
be reconciled with the decay rate found in these afore-mentioned calculations? 

In this work, we develop a different formalism, one that tracks the evolution of the 
quantum state directly in real time. It is the field-theoretic analog of the Wigner-Weisskopf 
(WW) [22] method in quantum mechanics and completely bypasses the need for the existence 
of asymptotic states. In this article we introduce the method, test it in Minkowski space time 
and apply it to two massless conformally coupled theories in de Sitter space time: Xcp 4 ; A</> 3 . 
The latter provides the simplest realization of a non-conformally invariant theory in which 
de Sitter expansion leads to vacuum instability in perturbation theory. 

Our major results are that indeed we find that the X(j) A vacuum is stable, as expected. 
The only time dependence the state acquires under time evolution is a phase corresponding 
to corrections to the vacuum energy. However, for non-conformally invariant interactions 
(e.g. cubic) the vacuum exhibits an instability to decay as a consequence of a time dependent 
wavefunction renormalization which diverges at late times, showing the infrared nature of 
this effect. This method allows for an unambiguous determination of the stability of the 
Bunch-Davies vacuum. Furthermore, it is a non-perturbative generalization of the transition 
amplitude one computes in perturbation theory, thus valid at late times. It is also closely 
related to the dynamical renormalization group method of resumming secular terms [11, 18], 
which has been used recently to deal with the secular growth of correlators in de Sitter space 
due to infrared effects[17](see also refs. [19-21]). 
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2 The Wigner-Weisskopf Method 



We start with a review of the Wigner-Weisskopf method in the case of Minkowski space time 
where there is time translational invariance. 

Consider a system whose Hamiltonian H is given as a soluble part Hq and a perturbation 
Hi: H = Hq + Hi. The time evolution of states in the interaction picture of Ho is given by 

i±Mt))i = Hi(t)Mt))l, (2-1) 
where the interaction Hamiltonian in the interaction picture is 

Hi(t) = e iHot Hie~ iHot (2.2) 

This has the formal solution 

|¥(t)>j = E/(i,to)|*(to)>j (2-3) 
where the time evolution operator in the interaction picture U(t,to) obeys 

: Hi(t)U(t,t ) . (2.4) 
J2Cn(t)\n) (2.5) 

n 

where \n) form a complete set of orthonormal states; in the quantum field theory case these 
are many-particle Fock states. From eq.(2.1) one finds the exact equation of motion for the 
coefficients C n (t), namely 

Cn(t) =~iJ2 C m (t)(n\Hi(t)\m) . (2.6) 
in 

Although this equation is exact, it generates an infinite hierarchy of simultaneous equa- 
tions when the Hilbert space of states spanned by {\n}} is infinite dimensional. However, 
this hierarchy can be truncated by considering the transition between states connected by 
the interaction Hamiltonian at a given order in Hi. Thus consider the situation depicted in 
figure 1 where one state, \A), couples to a set of states which couple back only to \A) 

via Hi. 

Under these circumstances, we have 

CaH) = -iJ2(A\H!(t)\K) C K (t) (2.7) 

K 

C K (t) = -iC A (t)(K\Hi(t)\A) (2.8) 
where the sum over k is over all the intermediate states coupled to \A) via Hi. 



i-uitM 

Now we can expand 

!*(*)>/ = 
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\ A ) /> N\ \ A ) 

\k) \k) 

(k\Hj\A) x (AIHjIk) 
Figure 1. Transitions \A) O \k) in first order in Hi. 

Consider the initial value problem in which at time t = the state of the system \^(t = 
0)> = \A) i.e. 

C A (0) =1,C K = 0. (2.9) 
We can solve eq.(2.8) and then use the solution in eq.(2.7) to find 

C K (t) = -i fi^Hjit'^CAit^dt' (2.10) 
J o 

rt 



where 



C A {t) = - [ E(t,t')C A (t')dt' (2.11) 
Jo 

S(t,t') = ^(A|fl>(t)|/6)(«|fl-j(0|A) (2.12) 



This integro-differential equation with memory yields a non-perturbative solution for the 
time evolution of the amplitudes and probabilities. Inserting the solution for C A (t) into 
eq.(2.10) one obtains the time evolution of amplitudes C K (t) from which we can compute 
the time dependent probability to populate the state \k), \C K (t)\ 2 . This is the essence of the 
Weisskopf-Wigner[22] non-perturbative method ubiquitous in quantum optics [23]. 

The hermiticity of the interaction Hamiltonian Hj, together with the initial conditions 
in eqs.(2.9) yields the unitarity condition 

J2\C n (t)\ 2 = 1 . (2.13) 

n 

In general it is quite difficult to solve eq.(2.11) exactly, so that an approximation scheme 
must be developed. We do this below and then compare the approximate results to both 
exact results when they are available as well as other approximation schemes. 

2.1 Markovian approximation 

The time evolution of C A (t) determined by eq.(2.11) is slow in the sense that the time scale 
is determined by a weak coupling kernel S. This allows us to use a Markovian approximation 
in terms of a consistent expansion in derivatives of C A . Define 

W (t,O = / Z(t,t")dt" (2.14) 
Jo 
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so that 

E(t, = jW (t, t'), W (t, 0) = 0. (2.15) 
Integrating by parts in eq.(2.11) we obtain 

J* £(i, t') C A (t') dt' = W (t, t) C A (t) - J W (t, t') ^pC A (t') dt'. (2.16) 

The second term on the right hand side is formally of fourth order in Hj and we see how a 
systematic approximation scheme can be developed. Setting 

/■*' 

W 1 (t,t')= W (t,t")dt", Wi(t,0) = (2.17) 
Jo 

and integrating by parts again, we find 

W (t, t') ^C A (t') dt' = Wi(t, t) C A (t) + • • • (2.18) 

leading to 

I E(t,t')C A (t')dt' = W (t,t)C A (t)-W 1 (t,t)C A (t) + --- (2.19) 
Jo 

This process can be implemented systematically resulting in higher order differential 
equations. Up to leading order in this Markovian approximation the equation eq.(2.11) be- 
comes 

C A (t) [1 - Wi(t,t)) +W (t ) t)C A (t) = (2.20) 

with the result 

C A (t) = e-tfiWW, £(t) = * -iW (t,t) [1 + W!(t,t) + ■■■} (2.21) 

Note that in general £{t) is complex. In particular if the time integral in eq.(2.21) is secular 
in time, the real part of £ yields a time dependent phase while the imaginary part of £ 
determines the time dependent decay since the sum rule eq.(2.13) requires that the norm 
of the state either remains the same or diminishes. The non-secular contributions to the 
integral yield the overall asymptotic normalization of the state, namely the wave-function 
renormalization constant. The connection with the resummation scheme provided by the 
dynamical renormalization group will be analyzed in appendix. A. 

Since both Wo, Wi oc Hj, the leading order solution of the Markovian approximation 
is obtained by keeping £{t) = — z Wo (£,£). However, if W\(t,t) features secular terms that 
invalidate the perturbative expansion at late time then the full expression in eq.(2.21) must 
be used. We will see in section(5) an important case where this situation emerges. 

In the Markovian approximation the coefficients C K (t) become 

rt , 
C K (t) = -i / (K\H I (t')\A)e- i fo WW dt'. (2.22) 
Jo 
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The non-perturbative nature of this solution is manifest by comparing this result to the usual 
perturbative transition amplitude in perturbation theory 

Ma^k = -i f (iAHi{t')\A) dt'. (2.23) 
J o 

This clearly corresponds to keeping Ca{P) = 1 in eq.(2.10), and is therefore valid only at very 
early times. 

2.2 Exact Solutions 

In order to understand the domain of validity of the Markovian approximation developed in 
the previous subsection, we take the time in this subsection to compare its results to those 
obtained in a situation where time translation invariance obtains. This allows us to solve 
eq.(2.11) exactly via Laplace transform methods. Thus, suppose that Hi is time independent 
in the Schrodinger picture. We take the states \A) and {\k)} to be eigenstates of Hq, 

H \n) = E n \n) . (2.24) 

We will assume {n\Hj\n) = by redefining the non-interacting Hamiltonian Hq to include the 
diagonal matrix elements of the interaction; this amounts to diagonalizing the perturbation 
to first order in the interaction. 

Using this in eqs.(2.7, 2.8), we find 

C K (t) = -i(K\Hj\A) f e« E «- E ^' C A (t>) (2.25) 
J o 

and 

/oo 

dJp{(J) e^A— ')(*-*') (2.26) 
-oo 

where the spectral density p{oj') is given by 

p(J) = Y^mH^H^ - J) . (2.27) 

K 

Introducing the Laplace variable s and the Laplace transform of Cyi(i) as Ca(s), with the 
initial condition CA(t = 0) = 1, we find 



C A (s) 

with solution 



s + I duo' , , - , 

S + 1{U)' ~ Ea) 



(2.28) 



C A (t)= I —C A (s)e st (2.29) 

-ioo+e ji7rl 
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where the e — > + determines the Bromwich contour in the complex s-plane parallel to the 
imaginary axis to the right of all the singularities. Writing s = i(u — ie) we find 

Auit 



c A (t) 



duo 
2^~i 



to — ie — f°° dio' 

J — oo 



u!+ui'— Ea— ie 



(2.30) 



In the free case where p = 0, the pole is located at oj = ie — > 0, leading to a constant 
Ca- In perturbation theory there is a complex pole very near uj = which can be obtained 
directly by expanding the integral in the denominator near uj = 0. We find 



- , I -~-AE A -z A u> + i-± 

UJ + UJ — ti A — It 2 



where 



p(u>) 



duo' ~ -AEa -ZAOJ + i^ (2.31) 

A£ -"i.^'(& < 2 - 32 > 

r A = 27rp(£U) (2.33) 
za = V T dJ ./ ( ^, 2 (2.34) 

and T 3 stands for the principal part. The term AEa is recognized as the energy shift while 
is seen to be the decay rate as found from Fermi's golden rule. The long time limit of 
C A (t) is determined by this complex pole near the origin leading to the asymptotic behavior 



where 



Z A 



C A {t)~Z A e 
1 

I- ZA 



d 



1 + za ' ~" dE A 
is the wave function renormalization constant, and 



[E A + AE, 



AE A = Z A AE A 
= Za Fa 



(2.35) 
(2.36) 

(2.37) 
(2.38) 



Now let's compare these results to those found via the Markovian approximation. With 
T,(t,t') given by eq.(2.26), to leading order in Hi we find 



£(t) 



Z(t,t')dt' 



UJ' 



1 



-i(u'-EA)t 



(2.39) 



so that 



f\{t')dt' = t f°° duo' P(UJ,) 



oc 
oo 



(E A - w') 
p(uj') 



dUJ ' (TP 

(E A - u?y 



sin(u/ — EA)t 
(uj'-E A )t 

1 - cos [(uj 1 - E A )t 



(2.40) 
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Asymptotically as t — > oo, these integrals approach: 

sin(fa/ — Ea) t 



-00 (Ea - u>) 



1 



(E A - u') 



dJ 



p(J) 



1-cos [(fa/-.E^)t] 



(i?A-^) 2 

Using these results we find that in the late time limit, 

Wi(f,t) -> -*a, 

leading to the result 



t^OO 7Ttp(E A )+V / 
J—i 



(E A - u 



1 I £(t') dt' ->■ -zA£^ t ft- z A Z A 

'0 2 



(2.41) 
( 2 - 42 ) 
(2.43) 

(2.44) 



where A^,r^,z A are given by eqs. (2.32,2.37,2.33,2.38,2.34,2.36). Prom this we read off 

C A (t) = Z A e~ iAE A 1 1 (2.45) 



where we approximated e~ ZA Za ~ 1 — za Za = Za in perturbation theory. This is in complete 
agreement with the asymptotic result from the exact solution eq.(2.35). 

Introducing this result into eq.(2.25) for the coefficients C K we find to leading order 

\{k\H t \A)\" 



\Ck(oo)\' 



(Ef - E K y + 



L A 



, E% = E A + AE A . 



(2.46) 



This can be interpreted as follows. If \A) is an unstable state, the states \k) with E K ~ Ea 
i.e. those nearly resonant with the state \A), are "populated" with an amplitude oc 1/Ta 
within a band of width Ta centered at E^. Furthermore 



K 



OO 



p(u>) 



{E* - fa/) 2 + 



1 



(2.47) 



where we used eq.(2.33). In this situation, unitarity entails a probability flow towards excited 
states. 

When 1^4 = 0, i.e. \A) is stable, we find 

\C A (t)\ 2 = Z\~(l-2z A ). (2.48) 



Using eq.(2.35) with Ta = in eq.(2.25), we see that 

„ N#/|^}| 2 



\c K (t)\' 



m-E K ) 



1 - cos [(Ea - E K )t 



(2.49) 
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In the asymptotic long time limit we find 

]T|C K (oo)| 2 ^2P f° dJ— P^l— = 2z A . (2.50) 

with |CU(oo)| 2 + ^Pjk^a |Ck(oo)| 2 = 1. This result clearly exhibits how the unitarity result 
in eq.(2.13) is upheld in the non-decaying state case. Most of the probability remains in the 
initial state, although a perturbatively small amount of it, « 2z "flows" to the excited states. 

Our conclusion from this comparison of results is that the Markovian approximation is 
trustworthy at late times, which makes it the perfect tool to study vacuum stability. We 
can also compare its results to those of a different approximation scheme, the dynamical 
renormalization group [11, 18], which resums secular terms to extract asymptotic late time 
behavior. In the interests of clarity we relegate this discussion to the appendix, but the result 
is that the Markovian approximation also agrees with DRG results. 

We now turn our attention to the field theoretic implementation of these ideas. 



3 The Quantum Field Theoretic Wigner-Weisskopf Method in Minkowski 
Space 

In this section we apply the ideas of section 2 to Minkowski space scalar field theories. We 
consider two separate situations. First, we take a massless scalar field with a quartic potential 
and compute the late time behavior of the vacuum persistence amplitude. We will then treat 
the case of a scalar cp coupled to another (lighter) scalar x via a cubic coupling ipx 2 - This 
will allow us to compute both the vacuum persistence amplitude as well as the single particle 
amplitude, from which a decay rate can then be extracted, following the procedures described 
above. 

3.1 Vacuum Amplitude in ^ 4 

We take the normal ordered interaction Hamiltonian to be 



H T (t) = A J d 3 x : cj) A (x,t) : , (3.1) 



where the interaction picture field is quantized as usual in a volume V 

<f)(x, t) = -1= ^= [ a k e~ i{kt - ls) + at e m-k-x)^ _ ( 3 2 ) 

k 

The normal ordering cancels tadpole diagrams and makes it so that to leading order in 
Hj, the vacuum is connected only to four particle states \k) = |l£ ; 1^; 1^ ; lg 4 ): 

(hv iffsi ifei hjmm = A i-f 5^ ( £ h) ; (3.3) 

v i=i ^{kik 2 k z ki)^ 
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we have subtracted the zero point energy of the vacuum so that Hq\0) = 0. This matrix 
element is depicted by the left hand side in figure 2. The kernel E(t,t') of eq.(2.26) is given 
by 



V 4 \6hk2hk4, 

ki k^ 

(3.4) 

where 

J (2tt) 3 J (2tt) 3 7 (2vr) 3 lefciifcafcalfci + fe + fcal 

This describes the i/iree /oop diagram on the right hand side of figure 2. 





(«|ff/|o> (o|fr/|«><«|ff 7 |o> 

Figure 2. y> 4 theory: left: transition matrix element (k\Hi\0) with the four particle state \k) 



Prom the results of the previous section we find 

C (t) = e~ zo e~ iAEot , (3.6) 

where 

AE = -\*V [fl^%— f j - ^ -r- 1 - ^ (3.7) 

y f = t ( 27r ) 16 (fciAfctelto + fc 2 + fcalj (£ 3 =1 h + \ki + h + h\j 

=» , , ,' , tat-; : j I tt, (3-8) 



=1 ( 27r ) 3 16 (fciM&lfci + fc 2 + h\) fe 3 =1 h + \h + k 2 + 

We also find that the vacuum is stable 

T = 2irp(u' = 0) = (3.9) 

as expected. 

AEq is recognized as the leading order C(A 2 ) correction to the vacuum energy in the 
normal ordered theory. The volume factor just reflects the extensivity of the vacuum energy. 

As we will compare the results obtained in Minkowski space-time with the case of de Sitter 
cosmology in the next section, it will prove illuminating to understand eqs.(3.7, 3.9, 3.8) in 
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more detail. According to the adiabatic theorem the exact ground state |0) is obtained from 
the non-interacting vacuum |0) by adiabatically switching on the perturbation Hj — > e _e '*' Hj, 
with e — > + 

|0> = [/ £ (0,-oo)|0> (3.10) 

where U e (t,to) is the time evolution in operator in the interaction picture with adiabatic 
switching-on of the interaction Hamiltonian. To 0(A) we find 

-o 



\0) = \0) -if £|K)(4H>(i)|0)e _6 l*l<ft (3.11) 

J — oo 



where \k) are the four particle states and (K\Hj(t)\0) is given by eq.(3.3). With Hj(t) 
e lHot Hje~ tHot the time integral can be done straightforwardly leading to 



1° (n\Hj(t)\0) e-l'l dt = -iX S® ( V h) - 



1 

,k 3 k 4 )^ + h + k 3 + k 4 + ie) 

(3.12) 

It is clear that the transition amplitude (n\U e (0, — oo)|0) has nothing to do with the decay 
of the vacuum state, which is stable in Minkowski spacetime, but rather with the "dressing" of 
the vacuum by virtual excitations; the true vacuum state is a superposition of many particle 
Fock states. Comparing eq.(3.12) with eq.(3.8) it is clear that 

K 

is given by eq.(3.8) and is a manifestation of unitarity. Furthermore the exact energy of the 
adiabatically constructed ground state is given by 

E=Eo+AEll= mjjhm (3 . 14 , 

<0|0) 

which up to 0(X 2 ) is given by eq.(3.7). This example clarifies that the Weisskopf-Wigner 
method provides a non-perturbative alternative to the adiabatic theorem of Gell-Mann and 
Low [24]. 

3.2 Vacuum amplitude and 1-particle decay in ipx 2 

The cubic coupling between two fields gives us the opportunity to compute a non-trivial 
decay rate using the WW techniques. We first compute the vacuum amplitude as we did in 
the quartic case, though in this situation instead of dealing with a three-loop diagram, the 
corrections come in at two loops. 

Thus, consider a massive field ip coupled to a massless field x with the interaction Hamil- 
tonian being 

Hi = X J d 3 x Lp(x) x 2 (x). (3.15) 
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The leading order vacuum to vacuum matrix elements are shown in figure 3 and we find 

POO 

£ (M')= / du 1 Po {J) e-^'^-^ (3.16) 



with 

Po (oj') = VX 2 
We also have that 



d 3 p f d 3 k 5{u' -E p -k-\k+p\) 



(2tt) 3 J (2vr) 3 2E p 2k2\k+p\ 



P 2 + 777,2 



Wo (t,t)= doj'^P-i 



00' 



(3.17) 



(3.18) 



Taking the long time limit of the leading order contribution, we find, using eqs.(2.41,2.42): 



i I £(t')dt' 
'o 



-iAEq t — Zq 



where 



(3.19) 
(3.20) 

note that AEq and zq are quadratically and linearly divergent in the ultraviolet respectively. 
This observation will be important when we analyze the vacuum to vacuum amplitude in de 
Sitter space-time. The late time behavior of the vacuum persistence amplitude is found to be 



-V I du —, z =V 



00 

-OO 



C (t) = Z e 



-iAE t 



Z = e 



1 - z - 



(3.21) 




k|#/|0) 




Figure 3. (p\ 2 vacuum to vacuum amplitude. Thin lines correspond to (p, thick lines to \- 



The diagrams that are summed by this result can be clearly interpreted: they are those 
that feature the leading dependence on the volume V at each order in A 2 . For example at 
order A 4 there are two types of contributions, two disconnected diagrams of the type depicted 
in figure 3 with a contribution oc (\ 2 V) 2 , and another disconnected diagram in which one 
of the internal lines features a one-loop self energy correction, with a contribution oc A 4 V 
which is subleading in terms of the volume factors and represents the C(A 4 ) correction to the 
vacuum energy. 

Next we turn our attention to the behavior of the single particle amplitudes C-^vit). To 
leading order in Hi, the one <p state is connected to states with two-x quanta. There are two 
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contributions to £(i,i') in this case: a disconnected one, describing a change in the vacuum 
state as well as a connected one which incorporates the contribution of the decay process to 
the self-energy of the (p field (figure 4) . 





Figure 4. ip \ 2 one <p particle amplitude. Thin lines correspond to if, thick lines to \- 



The two-x states for each diagram are different. For the vacuum diagram, we have \k) 



E(i,t / ) = S (i,O + Si(p,i,O 
where T<o(t,t') is given by eqs.(3.16,3.17) and 

/oo 
-oo 

with 



pi(p,u') = A 2 



d 3 k 5{oj' -k-\p-k\) 



(2tt) 3 2E p 2k2\p- k\ 
Following the same steps as in the previous examples we find 



i l £(t')dt' 
'o 



AE + A£i(p)J t--£t-[zo + zx(p)] 



where AEq and zq are given by eq.(3.20) and 



J-oo & p — 10 J- 



(Bp ~ " ) 2 



11%, It , Vi , 1* ) with k\ + &2 + ^3 = 0, while the connected diagram requires Ik) = 1 1 — , l x 

P ki k2 kz 
with p = k\ + k2 ■ 

It is straightforward to find in this case 



fci fca' 



(3.22) 



(3.23) 



(3.24) 



(3.25) 



(3.26) 
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The decay rate for <p — > 2x is found to be 



5(E p -k-\p-k\) 



leading to 

dj(t) = 2bZi(p) e - iA ^* e ~ iA ^i(p)* e -^* (3.28) 

The interpretation of this result is clear: AEq + AEi(p) is the total energy shift of the 
one particle state with respect to the unperturbed vacuum state, whereas AE\{p) is the 
energy shift with respect to the "dressed" vacuum state. Similarly, the total wave function 
renormalization ZqZ\(p) is the overlap between the bare single particle state obtained from 
the bare vacuum by applying a creation operator and the full "dressed" single particle state. 

The phase Ai?o t is the vacuum energy and is common to all many particle states. It can 
be renormalized away by introducing a counterterm Hamiltonian proportional to the identity 
operator 

H ct = -AE . (3.29) 

Then each many particle state features a term in the equation of motion for their amplitudes 
C n given by i AEq C n , which exactly cancels the phase from the vacuum energy (disconnected 
diagram). Furthermore, renormalizing the vacuum state and defining single particle states 
with respect to the full "dressed" vacuum state, can be achieved by defining the renormalized 
amplitudes 

C n (t) = C n (t)Z^ (3.30) 

leading to 

C l ,(t) = Z 1 (p)e~ iAE ^ t e- r -? t (3.31) 

This renormalization procedure cancels the contribution of the disconnected diagram 
which multiplies all many particle amplitudes by Zq e - jAE o* This is familiar from renormal- 
ized perturbation theory; the vacuum diagrams cancel between numerators and denominators 
in expectation values or correlation functions, since these correspond to the renormalization 
of the vacuum state. However, these disconnected graphs are present in the time evolution 
of the quantum states themselves. The redefinition of particle states by dividing the "bare" 
amplitudes by the amplitude of finding the bare state in the "dressed" state is an important 
ingredient in the adiabatic hypothesis of Gell-Mann and Low [24]. 

3.2.1 Particle decay and particle production 

The decay of a tp particle of momentum p results in the production of two x particles with 
k± , k% . The final result for the amplitude of the state with a pair of % particles can be under- 
stood from a simple argument which will be useful in the discussion of particle production in 
de Sitter spacetime. From eq.(2.8) the amplitude of the two % particle state is given by 

C 2x (t) = -i(l^,l^|fl>|lJ) e-^-^NDf^^)^. (3.32) 

J 
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Taking the amplitude C^tt') ~ 1 during the decay time scale t' < l/r^, the integrand in 
(3.32) is ~ 1 for the resonant wave-vectors k + \p — k\ ~ E p and the integral yields a secular 
term that grows linearly with time, while for wave-vectors well outside the resonant band, 
the rapid phase averages the integral out. Since the amplitude C]v(t') ~ for t' » 
the resonant wave-vectors for which the phase nearly vanishes lead to 



c 2x {t » i/t v ) cx l/r 



(3.33) 



Thus "resonant pairs" are produced linearly in time within a time ~ l/r<^ during which the 
amplitude C^e is nearly constant. This argument confirms the general result eq.(2.46), and 
will be useful in the interpretation of particle production from decay in de Sitter spacetime. 

3.2.2 On the wave function renormalization of the vacuum state: 

In anticipation of a study of the time evolution of quantum state in de Sitter spacetime we 
analyze in detail the time evolution of the vacuum state in the case in which all fields ip, x 
are massless, because we want to use this study in Minkowski spacetime as a guideline for 
interpretation of results in de Sitter cosmology. In this case So(t, t') can be calculated exactly 
by carrying out the three body phase space momentum integrals in eq. (3. 16, 3. 17), and we find 



E (t,0 



iX 2 V 



(4vr) 4 (t - f 



(3.34) 



where e provides a convergence factor that regularizes the ultraviolet behavior of the integrals. 
From this it follows that 



i I £{t')dt' 
'o 



, AV 
*2(4tt) 4 

from which we can extract in the long time limit 

X 2 V 



t 



1 



-^r-- + r 



X 2 V 



2(4vr) 4 e 2: 



ZQ 



2(4vr) 



(3.35) 



(3.36) 



This clearly displays the quadratic (1/e 2 ) and linear (1/e) ultraviolet divergences of the energy 
shift and wave function renormalization respectively. 
We also find 



and 



W (t,t') 



2(4?r) 4 

A 2 V 
'2(4vr) 4 ( 



1 



1 



(t-f- ief (t - ief 



1 + i 



t 



-ZQ + 0(1/*), 



(3.37) 



(3.38) 



this last result confirming the asymptotic behavior of eq.(2.43). A similar calculation shows 
that the wave- function renormalization for the single particle state, zi(p), is ultraviolet finite. 

It should be clear from the above discussion, and the comparison with the exact solution 
using the Laplace transform, that the vacuum amplitude Co(t) provides a resummation of the 
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vacuum diagrams of the type depicted in figures (2,3), whereas the amplitude for the single 
particle state C\(t) provides a Dyson- like resummation of the self-energy diagrams depicted 
in figure 4. 

This discussion and the comparison to known results in Minkowski space time explicitly 
clarifies that the WW method provides a resummation of the perturbative series akin to the 
Dyson resummation of self-energies insertion but directly in real time. For the vacuum case 
it sums (to lowest order) the two (^> 3 ) or three (^ 4 ) loop diagrams for the vacuum energy, and 
for single particle states it resums (to lowest order) the one- loop self energy diagrams (for cp 3 ) 
along with the vacuum diagrams. The loop corrections are given by the usual momentum 
integrals. 

Although the 3 theory is unbounded below and therefore has an intrinsic instability, 
such instability is not manifest either in the vacuum wave function renormalization or in 
single particle states for the self-interacting case (x = 4>) since kinematics prevents 1 — > 2 
processes. 

We will resort to these results when we try to interpret the various contributions to the 
evolution of the vacuum state in de Sitter spacetime. 

4 The Quantum Field Theoretic Wigner-Weisskopf Method in de Sitter 



The previous section showed quite clearly that the solution to the WW equations given by the 
Markovian approximation reproduces the results of standard flat space perturbation theory 
and and provides non-perturbative control of the late time behavior of the vacuum and many 
particle amplitudes C n (t). Having established the reliability of the non-perturbative result 
in Minkowski spacetime, we now turn to studying the time evolution of states in de Sitter 
spacetimes. 

We will consider the same situations as described in the previous section, i.e. first that of 
a quartically self-coupled field and then a cubically coupled one. The new features to be added 
to the discussion are first, that we will consider only massless, conformally coupled fields (i.e. 
£ = 1/6), and second, in the cubic case we allow for the possibility that \ = ip. In Minkowski 
spacetime this last condition would preclude any decay phenomena from occurring, but as 
we shall see below, the lack of a global time-like Killing vector in de Sitter spacetime allows 
energy non-conserving processes to occur. 

We consider a spatially flat Friedmann-Robertson- Walker (FRW) cosmological spacetime 
with scale factor a(t). In comoving coordinates, the action is given by 



Space 




(4.1) 



with 




(4.2) 
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being the Ricci scalar, £ = corresponds to minimal coupling while £ = 1/6 corresponds to 
conformal coupling. 

Specializing now to the de Sitter case with a(t) = e Ht , it is convenient to pass to conformal 
time r] = —e~ /H with drj = dt/a(t) and introduce a conformal rescaling of the fields 



a(t)<f>(x,t) = x{x,rj). 



(4.3) 



The action becomes (after discarding surface terms that will not change the equations of 
motion) 



s 



1 



d 3 x dr\ 



1 



(4-p) 



with primes denoting derivatives with respect to conformal time r\ and 



M 2 Jr,) = (M 2 + ZK )C 2 ( V ) 



where for de Sitter spacetime 



Cirj) = a(t( V )) 



1 



In this case the effective time dependent mass is given by 



M 2 J V ) 



M 2 / 1 



(4.4) 

(4.5) 
(4.6) 

(4.7) 



The Heisenberg equations of motion for the spatial Fourier modes of wavevector k of the fields 
in the non-interacting (g = 0) theory are given by 



k 



V 



W 2 1 



Xjfrl) = > 



where 



■ M 2 

I — + 12£ 



, 2 = 9 --(t 



Two linearly independent solutions are given by 

g v {k;rf) = i r v ~\ ^fj H® (k V ) 

f„(k; V ) = \ r+Vy/mjHQ-Hkri) = [gv(k; V )}* , 

(1 2) 

where Hu ' (z) are Hankel functions. Expanding the field operator in this basis 



V 

k 

The Bunch-Davies vacuum is defined such that 

a fe H0) = 0, 



(4.8) 
(4.9) 

(4.10) 
(4.11) 

(4.12) 
(4.13) 
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and the Fock space states are obtained in the usual manner, i.e. by applying creation operators 



at to the vacuum. 



In the Schroedinger picture the quantum states |^(??)) obey 

ijj-\y(r)))=H(r))\*( V )) (4.14) 

where in an expanding cosmology the Hamiltonian H (77) is generally a function of 77 (unlike 
in Minkowski space-time where it is a constant by energy conservation). Introducing the time 
evolution operator U(rj, 770) obeying 

i-jj-Ufa Vo) = H( V ) U( V , r/ ), ^(770,770) = 1, (4.15) 

the solution of the Schroedinger equation is \^f(ij)) = {7(77,770) \^(t]q)). Writing the Hamilto- 
nian as H(rj) = Hq(i]) + Hj(rj) with Hq(t]) the non-interacting Hamiltonian, and introducing 
the time evolution operator of the free theory Uq^^tjq) satisfying 

d d 
i-j^U (rj,ri ) = H (rj)U (r),r) ), i—ZJ^" (77,770) = -Uq (v,m) h o(v), Uo(rjo,Vo) = 1, 

(4.16) 

the interaction picture states are defined as 

|*(r?))/ = f/ - 1 (r?,77o)|^(r7)). (4.17) 

These obey 

^(V))i = -tfliWItffa))/, H T (ri) = f/ - 1 (7?,r / o)^(77)f/o(r?,%) (4.18) 
where the normal-ordered interaction Hamiltonian in the interaction picture is 

HM = { _ H ^- P) j ■■ \X& VW ■■ d 3 x (4.19) 

where x is the free field Heisenberg field operator in eq.(4.12). We now expand the state 
^(77))/ in Bunch-Davies Fock states \n) \^(r]))j = J2 n Cn(jl)\ n )- 

We can now implement the Weisskopf-Wigner program directly by replacing t — > 77 in 
the results of the previous sections, with the proviso that t — > 00 => 77 — > 0~, and the initial 
conditions on the states are given at some given initial time 770. Eventually we will take 
770 — > —00. Equations (2.10,2.11) now become 

C K {rj) = -i r(K\Hitf)\A)C A {rf)dr{ (4.20) 



Ca(v) = ~ ( V S(r7, 77') C A (rf) drf (4.21) 
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where 

= Y,( A \ H dvMm I (v')\A) (4.22) 

K 

and the intial conditions Ca(vo) = 1 ; C k ^a(vo) = 0- 

In the Markovian approximation to leading order in the interaction 



C A (n) = e -!l w ^'^^\ W (v',v') = f nv',v")dv" ■ (4.23) 



We now specialize to the case described at the beginning of this section, namely M 2 = 
0, £ = 1/6, v = 1/2, in which case canonical quantization of the scalar field in a finite 
comoving volume V leads to 



X(x, n) = -^L V -L \a r e~ ikri e** + at e ikr > e~ ilg ] . (4.24) 



r 

is 

4.1 De Sitter Vacuum Amplitude in 4 



A massless, conformally coupled quartically self-coupled scalar theory is conformally (or more 
precisely Weyl) invariant and to lowest order (i.e. neglecting higher order radiative correc- 
tions that break conformal invariance via renormalization effects), correlation functions are 
independent of the metric. 

The vacuum to vacuum amplitude is described once again by the diagram in figure (2) 
leading to 

s o (??,?/) = x v / n44 1gl ; ; ,r , r ■ r, ^ 



" ( 27r ) lQk^hlh +k 2 + k 3 \ 



1 _ e -i(T,Liki + \ki+k 2 +k 3 \)(r,'-r l0 ) 



Wotf, rf) = -i\ 2 V I f[ 44 U ^ =, A ^r. (4.26) 

J f = \ (2vr) 3 mhhhlh + k 2 + k 3 \ [ E-=i fci + \h + k 2 + fe 3 |] 

Integrating in 7/ and taking 77 — 770 — cxd we find 

C (7?) = Z e -i*Zo(v-vo) (4.27) 

where i?o ~ 1 ~~ z o an d AEq are the same as the Minkowski spacetime results of eqs. (3.7,3.8). 

As expected due to the Weyl invariance of the system, the vacuum is stable. There is 
only an n dependent phase and an overall wave function renormalization just as in Minkowski 
spacetime. However, the main point of this exercise is to clarify that a transition matrix 
element {k\Hj\0) 7^ is not a signal of any instability of the vacuum (or for that matter of 
particle states). Rather, it points to the "dressing" of the particular state via its coupling 
to a continuum of states and the ensuing renormalization. Only the time dependence of the 
amplitude can distinguish between a true decay of the state or a simple phase and this is 
exactly what the Wigner-Weisskopf approach focuses on. 
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4.2 Vacuum amplitude and 1-particle decay in i^ 3 

A more interesting situation appears when we consider an interaction that is not Weyl in- 
variant, such as a cubic self-interaction tp 3 . The reason for this is that as seen in eq.(4.4), the 
interaction vertex for x now depends on the scale factor, and hence on the conformal time r\. 

We follow the approach we took in the flat space case, first concentrating on the vacuum 
amplitude. In the Markovian approximation, the vacuum persistence amplitude satisfies 



(4.28) 



The matrix elements for So are depicted in figure 3 but with all the lines now corresponding 
to the same massless, conformally coupled field ip. We find 



£o(W) 



\ 2 V 



d 3 p f d 3 k e~ i ( p+k+ ^ + ^( v ~ vn> 



H 2 rjr]' J (2vr) 3 y (2vr) 3 2p2k2\k + p\ 



(4.29) 



The momentum integrals are similar to those appearing in the Minkowski case with massless 
fields and yield 



(4vr) 4 H 2 nr]' 



Tj — Tj 



It 



Here e — > + provides a ultraviolet convergence factor for the momentum integrals just as in 
Minkowski spacetime. We obtain 



W (77,77) 



X 2 V 



(4tt) 4 H 2 rf 



V 

2e 2 



In 



-V 



n ir 
+ 7 + 2 



where the dots stand for subleading terms as ry, e — > 0, r/o — > — oo. 



find 



where 



A (r?) 



C (v) - e 



X 2 V 



-iAo(Tj) e -zo(»?) 



2e 5 



2(4 7 r) 4 tf 2 (-7 ? )e 2 
A 2 V 



3 rj 2 



In 



2(4vr) 4 ff 2 r? 2 e 



1 + 



7r e 
3rj 



+ ... 



Keeping the leading order terms in the asymptotic limits n — > , 770 



-00 , e 



(4.31) 

0+ we 
(4.32) 

(4.33) 
(4.34) 



The 1/e 2 and 1/e ultraviolet divergences are the same as in Minkowski spacetime, when 
we match t — >■ 00 to rj — > 0~. However, there are some important differences with the flat 
space result. The wave function renormalization, featuring the 1/e divergence, is now time 
dependent and grows as rj — > CP. This is a consequence of the extra factors 1 / 77 in the 
interaction vertices. The phase Ao(^) also features a novel logarithmic ultraviolet divergence. 
The time evolution of the wave function renormalization function zq (77) is such as to force 
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the decay of the vacuum state and this decay is hastened as rj — > 0~, revealing the infrared 
nature of instability of the vacuum. 

Comparing the wave function renormalization given by eqn. (4.34) with the result in 
Minkowski space-time given by eqn. (3.36) it is clear that the instability has nothing to do 
with the intrinsic instability of the </> 3 theory from being unbounded below but is a direct 
manifestation of the cosmological expansion and breaking of conformal invariance, since the 
conformally invariant features an ultraviolet divergent but r\ independent wave-function 
renormalization constant. 

The leading order ultraviolet divergence of Zo(rj) acquires an interesting interpretation if 
we assume that the short distance/time ultraviolet cutoff e is constant in physical coordinates: 

constant. (4.35) 



(-H V ) 

Under this assumption, the leading order term in zq becomes 



, A A 2 V phys (rj) 
MV) = 2(47r)4? (4-36) 

where 

is the physical volume. This is a striking result: the vacuum wave function renormalization 
with this regularization is the same as in Minkowski spacetime, eq.(3.36), but in terms of 
the physical volume which grows as a function of time. Therefore with this choice of regular- 
ization, the decay of the vacuum state is a consequence of the fact that the spatial volume 
increases with time, hence the amplitude for finding the bare state in the dressed state falls 
off exponentially with the physical volume. 

This behavior of the vacuum state leads to some rather interesting consequences for 
various correlation functions that might be associated with physical observables, such as non- 
gaussianity. We do understand that inflaton fluctuations £ behave as massless minimally 
coupled fields with derivative and non-local interactions, and so our calculations will have 
limited utility in understanding what happens with £ correlators. However, we expect that 
the WW method can be brought to bear on this question and so we present here a calculation 
of the effects of the decay of the vacuum on the non-gaussian expectation value 

<0| :x 3 (x,?7) : |0>. 

Passing to the interaction picture, xO?j v) is given by eq.(4.24) and \0(rj))j = Yl n ^n{v)\ n ) 
where the coefficients C n {rf) are the solutions to the Weisskopf-Wigner equations (4.24,4.21) 
with A = 0. In an obvious notation : x 3 (x, r/) :~ (o^) 3 + (a) 3 + (a^) 2 a + a^a 2 and since to 
lowest order Hi connects the vacuum to the three particle state \n) = |1^; 1^; lj£ )> (see 
figure 3), the terms that contribute are of the form 
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The coefficients can be computed using 

(i^i^jh^m = :„::y.7i' (4.38) 

which leads to 



A (2,)3^)( fcl+fc2 + fc3 ) Aki+k2+k3)v 

fe2 ' * 3 ' ' Km ' Hr, l8V3 kl k 2 k^ 



CM = -ij-W^^^) f e^^W ^ *£. (4 . 39) 



Summing all contributions we finally find 



cm r e^+fe+fcow) Co(r? ') ^ +h . c . . 



-n J (2ir) 6 8kik 2 \h + k 2 \ Jr, n 

(4.40) 

The lowest perturbative order would replace Co — > 1 leading to an infrared logarithmic di- 
vergence as 77—7-0. However the Weisskopf-Wigner resummation replaces Co by the solution 
(4.32) thereby suggesting that vacuum decay severely damps out the non-gaussian correla- 
tions. 

While this is a very preliminary result, it does suggest that the resummation provided 
by the Weisskopf-Wigner method may lead to important modifications of non-gaussian cor- 
relations, a topic of much importance [19-21] that will be explored elsewhere. 

Now let's turn to the behavior of the one-particle states here. Recall that in the Minkowski 
case, the Markovian approximation to the WW solutions allowed us to extract a decay rate 
for the state, as well as an energy shift and a wavefunction renormalization. The situation in 
de Sitter is, as might be expected by now, both more complicated and far more interesting. 

For single particle states there are two contributions to £(77,7/) just as in the single 
particle case in Minkowski spacetime: £(77,77') = £o( 7 ?) 7 ?') + ^liPiViV')- The disconnected 
vacuum contribution Eo(f7, ?/) is given by eq.(4.30) whereas for the connected contribution 
we find 

A 2 f d 3 k e -*(k+|p-£|-p)('7-'7'- ie ) 



where e — > + is a convergence factor that regulates the ultraviolet behavior of the momentum 
integral. This integral can now be carried out leading to 

-i\ 2 

Si (p, 77,770 = (4.42) 

S2ir z p H A 77 77' (77 — 77' — le) 



from which we obtain the connected contribution to Wq- 

, „ , , i A 2 

' no 



W ,i(p, r?, 77) = £ E!(p, 77, r/)drf = 32 J* H2 ^ j - i\ - In (-?) - In ^ 

(4^43) 

The ln(e) reflects a logarithmic ultraviolet divergence, also seen in the one-loop self-energy, 
which can be renormalized by a mass renormalization. Including a mass counterterm in the 
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Lagrangian 5M 2 : (f) 2 : /2, yields an extra term in the interaction Hamiltonian. Upon passing 
to conformal time and rescaling the fields by the scale factor we have 

Rf {r]) = 2Ww I * x SM * :x2 '" (444) 

This yields the matrix element 

-z(l p H^(r ? )|V) = -,^- f . (4.45) 
Taking this counterterm to first order and choosing 

5M 2 = * In (1) (4.46) 



16tt 2 v/i> 

where ji is a renormalization scale, results in renormalizing Wo i by replacing e — > ji in 
eq.(4.43). 

After mass renormalization we find the leading asymptotic behavior as rj — > 0~, r/o — > —oo 

C v (r?) = (70(7?) C 1? (ry) , C v (r?) = e~ iAl ^ e~^) ( 4 .47) 

where 60(77) is given by the vacuum (disconnected) contribution eq.(4.32) with eqs. (4.33,4.34) 
and Ci-(r?) is the one particle connected amplitude with 

A 2 In (=) 

= (448) 

The decay law e" 7 ^) with eq.(4.49) can be understood without invoking the Markovian 
approximation as follows: Si(p, 77,7/) given by eq.(4.42) can be written as 

Si(p,W)= o 9 2~ t ff2 > V T^—i\ + io X m 2 ^-V)- (4-50) 
S2ir z p H z rj rj' (rj — rj') S2irpH z rj z 

The imaginary part leads to an overall //-dependent phase for Ci, while the real part leads to 

+ aA A L 2 Ci(r?) = Ci(t/) = (phase) x e~^) (4 . 51) 
o47r pH z rj z 

The extra factor 1/2 arises because the argument of <5(r/ — rj') vanishes at the upper limit of 
integration. 1 Thus the asymptotic decay law is valid beyond the Markovian approximation. 
An as already pointed out in ref.[10, 11] the decay of single particle states into pairs of their 
own quanta is a consequence of the lack of kinematic thresholds in the expanding cosmology. 

The single particle decay 1 — > 2 in the self-interacting massless (f) s theory in de Sitter 
space time is again in striking contrast with the situation in Minkowski where kinematics 
forbids this process. Therefore this single particle instability is again a consequence of the 
cosmological expansion and the manifest lack of kinematic thresholds and not in any manner 
associated with the intrinsic instability of 3 mentioned above. 



lr rhis can be seen by writing the integral as f° Q(r] — rj')5(rj — rj') Ci(r/) drj' with 0(0) = 1/2. 
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4.3 Pair Production 

Just as in Minkowski spacetime, the decay of the single particle state entails the production 
of correlated pairs of particles with momenta k, p — k. The amplitude for pair production 
is obtained from the general renormalized form eq.(B.5) (see Appendix B for details) where 
\A) = \h 



•-pn 



C(k,q,rj) 



i A (2vr) 3 6W(k + q-p) f n 



H 



8V 3 pkq 



J,(k+q-p)r)' n (I) d rf 
rf 



(4.52) 



The rf integral in eq.(4.52) cannot be done exactly; however, for weak coupling 2 , X 2 /H 2 <C 1 
we can provide a reliable estimate of the amplitude. 
In the integrand in eq.(4.52) 



din) e F55 

- (phase) x — — , g 



X 2 



-T] 



(-v) 



64-irpH 2 ' 



(4.53) 



This function is oc l/( — 77) for — r/ > g, reaches a maximum at — 77 ~ g and then falls off sharply 
to zero for —r/<g. This means that the integrand is dominated by the region — r\ > g in 
which we can approximate C\ ~ 1. 

Approximating C\ p ~ 1 for —r/ > g in eq.(4.52) and assuming that all wave-vectors are 
deep inside the horizon at the initial time 7/0, namely —(k + q — p)r/o 3> 1 we find 



C(k,q;ri) ~ i 



X(2Tr) 3 5^(k + q-p) 



H 



8V 3 pkq 



■Ci[-Kr/]+iSi[-K7/] 



7T 



K = k + q-p (4.54) 



where Ci, Si are the cosine and sine integrals respectively. It is clear that the amplitude for 
pair production is largest when K becomes superhorizon, namely —Kr/ 1, in which case 
taking —r/ ~ g, we find 



C(k,q;r]) ~ i 



A(2vr) 3 5^ 3 \k + q-p) 



H 



8V 3 pkq 



In 



1 

YKg 



7T 



Kg < 1. 



(4.55) 



For modes that remain well inside the horizon during this time scale, the oscillatory phase 
leads to strong dephasing and cancellation of the integral since Ci[— Krf\ — > 0, Si[— Krj\ — > it/2 
for — Kr] S> 1. 

We see then that the decay of single particle states leads to the production of particles 
with superhorizon wavelengths. The total probability can be estimated by summing over 
the superhorizon wave-vectors that yield the largest contribution, namely those for which 



2 In the X(f> 3 theory A has dimensions of mass and the perturbative expansion is in terms of the dimensionless 
ratio \/H. 
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K ^ Vs- Since in perturbation theory gp = X 2 /QAttH 2 <C 1, it follows that these wave- 
vectors correspond to k < 1/g. An approximate estimate for the total probability of pair 
production is obtained by neglecting the logarithm in eq.(4.55), leading to 



^\C(k,q ]V ~i/g)\ 2 



A 2 



I 



9 



p + k — \p — k\ dk 



(4.56) 



32vrV# 2 



k,q 



for gp <C 1. The integral in eq.(4.56) ~ 2p/g leading to the approximate result 



^IC^g;^!^)! 2 



2 



(4.57) 



1 



7T 



k,q 



This result is a confirmation of unitarity. 

It is important to establish a difference with the case of single particle decay in Minkowski 
spacetime analyzed in section (3.2). Particle decay in Minkowski spacetime leads to pairs 
with a non-perturbatively large amplitude oc 1/T but in a narrow band of wave- vectors with 
a width oc V. Unitarity is fulfilled by populating a narrow band of nearly resonant wave- 
vectors with large amplitudes. In de Sitter space-time the situation is very different; the 
massless particle of momentum p decays into particles with wave- vectors k,p — k which are 
non-resonant. The production amplitude begins to grow logarithmically when k crosses the 
horizon with p either sub or super horizon. The analysis above demonstrates that particles 
with superhorizon wavelengths are produced with perturbatively small amplitudes, but at long 
time more and more particles become superhorizon and these states all become populated 
up to wave- vectors of order 1/g. The sum total of probabilities is very nearly one, up to the 
logarithmic corrections that were neglected in the calculation above. 

The detailed process of particle production actually depends on the time scale in a subtle 
manner. If the wave-vector p is deep inside the Hubble radius during some early time scale, 
namely — pn 3> 1, the condition for the phase to be slowly varying and thus for the integral 
to grow logarithmically, \{k + \p — k\ — p)w\ <C 1, implies p 3> k and \kn\ <C 1, since for this 
configuration 



This means that the decay process corresponds to the emission of soft superhorizon quanta 
and a recoil of the particle, just like bremsstrahlung as found in ref.fll] for the decay of sub- 
horizon modes. After p becomes superhorizon, pair production becomes more efficient by 
the decay into nearly back-to-back superhorizon entangled pairs of momenta p <C k < 1/g. 
Therefore particle production hastens when p crosses the horizon and this decay process leads 
to correlated pairs of superhorizon particles. 

From the definition of g in (4.53) it follows that pg = \ 2 /QAirH 2 <C 1 for weak coupling. 
Therefore, in the weak coupling regime the wave-vector p exits the Hubble radius well before 
the time scale — n ~ g at which C\ p begins to decrease and production of superhorizon pairs 
occurs during a large interval in time during which the amplitude C\ ~ 1. 
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We conclude with the observation that at long times decay of single particles in de Sitter 
spacetime leads to an entangled state of superhorizon pairs up to momentum ~ 1/g, 

1/3 _ 

\lp)^Y,C(k,p-k;n~l/g) \1 R ,1^) (4.58) 
k 

where the pair amplitude is given by eq.(4.55). Since p 1/g and the integral in eq.(4.56) 
is dominated by momenta near the upper limit ~ 1/g, it follows that the processes with the 
largest amplitudes correspond to the emission of nearly back-to-back superhorizon pairs with 
momenta k^> p. 

5 Non-perturbative self-consistent screening? 

The above results were obtained under the assumption that W\{t,t) in eq.(2.21) is perturba- 
tively small and can be neglected. While it is formally true that W\{t,t) oc A 2 , the results 
for the vacuum eqs. (4.32-4.34) might give us pause as to the validity of this argument. The 
volume, cutoff and 77 dependence which becomes singular at long time 77 — > 0~ all conspire 
against the validity of perturbation theory. However, if zq does not depend singularly on n 
in the long time limit and reaches a constant, as is the case in Minkowski space time then 
this is simply interpreted as a (large) renormalization and "dressing" of the bare state, which 
is renormalized by wave function renormalization as is the case in the adiabatic hypothesis. 
However, the strong n dependence of zo(r]) in (4.34) clearly indicates that the fully "dressed" 
vacuum state becomes orthogonal to the bare vacuum state at long time. Thus unlike the 
case of Minkowski space time where W\ approaches a constant and its contribution can be 
simply lumped together in with the constant (albeit cutoff and volume dependent) wave func- 
tion renormalization, in de Sitter space time Wi (77,77) grows singularly as 77 — > 0~. Thus it 
is possible that the denominator in eq.(4.32) compensates for the growth of the numerator 
slowing down the decay of the vacuum state. To assess this possibility we now calculate W\ 
in eq.(2.21): 

where again we have neglected subleading terms. As in Appendix B, we allow for a countert- 
erm H c t(r]) in the Lagrangian to cancel the (infinite) vacuum phase, so that the equation for 
Co now becomes 

C (7?)[l - ^1(77,77)] + [iH ct ( V ) + Wo(v,v))C {v) = (5.2) 

with Wq (77,77) given by eq.(4.31). Requiring that H c t cancel the divergent imaginary part of 
Wo(r], 77), we write 



G 



_ ln(-77/e) 
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defining the variable $(r/). Remarkably, a suitable choice of $(77) leads also to a cancelation 
of the real part of Wq , namely the contribution that leads to vacuum decay. With the choice 
of counterterm eq.(5.3) we find 



iH ct + W (rj,rj) _ 



where 



77 7T 

a ~ ~2~e ~ 2 



(5.4) 

(5.5) 
(5.6) 



In the limit G/(—n) s 3> 1 the ratio in eq.(5.4) becomes 



+ 



$(r7)/3(r7)-«(r?)(i-^) 



a 2 (7?) + /3 2 (77) 



(5.7) 



The real part of the above expression leads to the time dependent vacuum wave-function 
renormalization and the decay of the vacuum state. However, it is noteworthy that the 
particular choice 



(5.8) 



cancels the real part leading to 



iH ct + W (n, 77) 
1 - ^1(77,77) 



2'/ 



In [=2 

e e 



(5.9) 



namely a purely imaginary contribution. With this choice of $ the solution of eq.(5.2) becomes 



Co(r?) 



0(7?) 



eln =? 



+ 0(ln(ln(-77))) 



(5.10) 



We interpret this result as a self- consistent screening of the vacuum evolution. The result 
in eq.(5.4) is strongly reminiscent of the behavior of a running coupling towards a fixed point; 
in d = 4 — e (Euclidean) dimensions the dimensionless renormalized dimensionless coupling 
of a A<^> 4 theory approaches an infrared fixed point, namely 



1 + 



Are 



* 7 



(5.11) 



where I is a constant and k an infrared scale. This suggests our conjecture: that non- 
perturbative self-consistent effects may lead to a fixed point in Hilbert space in which the 



-27- 



vacuum state is stable, evolving in time only with a phase. In this analysis we only focused on 
the terms that diverge as rj — > 0~ , e — > and neglected terms that remain finite in these limits 
which result in a volume dependent but finite contribution to the phase and wave function 
renormalization which remain perturbatively small in these limits. 

At this point we can only conjecture this behavior based on the above analysis, and 
clearly much more needs to be understood before concluding that this effect is generic to 
all theories. Furthermore, it could well be that higher order effects that may have been 
missed could change this picture. An important aspect that will be discussed elsewhere is the 
manifestation of unitarity from the contribution of states that are connected to the vacuum 
by the interaction. Clearly, this phenomenon is noteworthy of further investigation and its 
impact on the asymptotic behavior of correlation functions require further and deeper study. 

6 Conclusions and further questions 

In this article we introduced a non-perturbative formulation to study the time evolution of 
quantum states in cosmological space time. The method is a quantum field theoretical gen- 
eralization of the Weisskopf-Wigner[22] theory of atomic line widths which leads to a resum- 
mation of the perturbative series for the real time dynamics. In a Markovian approximation 
to the equations for the amplitudes, we established the correspondence between this method 
and Dyson-type resummations in Minkowski spacetime as well as with the dynamical renor- 
malization group approach to real time dynamics [18]. In Minkowski spacetime this method 
allowed us to establish the asymptotic long time limit of the amplitudes of quantum states 
thereby revealing how unitarity is manifest both in the case of non-decaying vacuum and 
decaying single particle states. Having established the reliability of the method in Minkowski 
space-time, we generalized it to study the time evolution of quantum states in cosmological 
space time. 

We then began a program to study long time and infrared divergences by implementing 
the non-perturbative method in two prototype quantum field theories in de Sitter space time: 
conformally coupled massless (j) 4 , <p 3 . These are the simplest scalar field theories that still 
possess many of the relevant features that are of current intense interest. The t/> 4 theory 
is conformally invariant and the asymptotic dynamics of quantum states is similar to that 
of Minkowski spacetime: the vacuum is "dressed" but stable with an amplitude having a 
divergent but time independent wave function renormalization and a phase which is also 
divergent and reflects vacuum fluctuations. 

The (p 3 case is more interesting; even when conformally coupled to gravity and massless, 
it is not a conformal invariant theory and particle production takes place. We found that the 
vacuum amplitude features not only a divergent phase, but also a wave function renormaliza- 
tion which depends singularly on (conformal) time and is proportional to the volume. This 
behavior entails an instability of the vacuum state which is enhanced by the volume factor. 
We find that in a particular regularization scheme, this instability can be simply understood 
from the vacuum wave function renormalization in Minkowski spacetime, but in terms of the 
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expanding physical volume thereby arguing that the decay of the vacuum is a consequence 
(in this renormalization scheme) of both the wave function renormalization of the vacuum 
state being proportional to the spatial volume as well as the cosmological expansion of the 
physical volume. The potential impact of vacuum decay on a simple non-gaussian correlation 
function was discussed. We showed that the vacuum diagrams in the time evolution of excited 
quantum states can be canceled systematically in perturbation theory by renormalizing the 
amplitudes with a time dependent vacuum wave function renormalization. In this theory we 
also showed that single particle states decay into two quanta of the same field, as a conse- 
quence of the lack of kinematic thresholds, resulting in particle production. This process is 
hastened when the momentum of the decaying particle becomes superhorizon leading to an 
entangled state of superhorizon pairs. We also analyzed how the amplitude of quantum states 
of these produced pairs manifestly satisfies unitarity. 

The singular time dependence of the vacuum wave function renormalization, combined 
with the volume and cutoff dependence all suggest a breakdown of the leading Markovian ap- 
proximation. We considered the general expression eq.(2.21) for the Markovian approximation 
and found the remarkable result that a particular choice of counterterm in the Hamiltonian 
leads to a complete cancelation of the secular and singular contributions to the vacuum wave 
function renormalization at long time and the asymptotic stability of the vacuum state. We 
conjecture that this is a self-consistent "screening" mechanism akin to the approach to a fixed 
point of a running coupling, in this case describing a dynamical fixed point in Hilbert space 
in the evolution of the state. 

This body of results lead us to ask further questions that will be addressed in forthcoming 
work: 

• How to systematically calculate correlation functions in the quantum states whose am- 
plitudes are solutions of the generalized Weisskopf-Wigner equations. This translation 
between the dynamics of states and correlation functions will allow a direct comparison 
with the dynamical renormalization group resummation program implemented at the 
level of correlation functions in refs.[ll, 17]. 

• How to consistenly go beyond the Markovian approximation studied here, and in partic- 
ular, how to understand whether going beyond the Markovian approximation introduces 
important corrections to the results obtained here. 

• In minimally coupled massless theories, there are additional infrared divergences. How 
do these manifest themselves in the resummed amplitudes for the quantum states? Is 
there an infrared fixed point that emerges from the resummation program? 

• Ref. [25] provides an intriguing, all orders result stating that the vacuum of a massive 
scalar field in de Sitter space time is both stable and free of infrared singularities, and 
that correlation functions may be computed in a Euclidean compactification and ana- 
lytically continued to Lorentzian signature without instabilities or infrared divergences. 
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These results seemingly contradict those of ref . [9] . Clearly it is important to understand 
whether the asymptotic screening mechanism found in our work can explain the results 
of ref. [25] within the real time dynamics of states, or at least provide a bridge towards 
a translation between the Euclidean and real time formulation at asymptotically large 
time. 

Work on these and other questions is in progress and will be reported elsewhere. 
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A The DRG and the Markovian approximation 

The DRG approach also offers a systematic method to study the evolution of quantum 
states [18] and provides an alternative to the adiabatic switching-on procedure to generate 
exact eigenstates of the interacting Hamiltonian from the eigenstates of the non-interacting 
system. Along the way this DRG formulation extracts the energy shifts, wave function renor- 
malization constants and offers a clear description of the decay dynamics of quantum states 
and its description in terms of Fermi's Golden rule. 

Consider the time evolution of a quantum state in an interacting theory with Hamiltonian 

H = H + XH I , (A.l) 

namely 

= e -mt-to) = e -iH t u{tjto) e iHoto mo)) f (A . 2) 

where we have introduced the unitary time evolution operator in the interaction picture 

ft rt ft' 

U{t,t ) = l- i X / H^t') dt' - A 2 / / iTj(t') fl>(t") dt' dt" + C(A 3 ) , U(t ,t ) = l, 

J Iq J to J to 

(A.3) 

with 

#/(t) = e iHot Hi e~ iHot . (A.4) 
It is convenient to pass to the interaction picture in which the states are given by 

W))i = e iHot \m) , (A.5) 
and their time evolution is given by 

m)) i = u(t,t )m )) i . (a.6) 
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Consider the adiabatic time evolution operator U e (t,to) obtained by replacing 

Hl (t) -> e- e l'l H T (t) , (A.7) 

in the time evolution operator U(t,to) in eq.(A.3). The adiabatic, or Gell-Mann-Low theorem 
asserts that the states 

|tf„> = C/ e (0,-oo) |n) , (A.8) 

constructed out of the eigenstates |n) of the non-interacting Hamiltonian Hq with energy E n , 
are exact eigenstates of the total interacting Hamiltonian H, and the exact eigenvalues are 
then given by [24] 

The dynamical renormalization group provides an illuminating alternative to this procedure. 
To make this approach more clear and to establish contact with familiar results, we now 
consider the case in which \ip(to))i is an eigenstate of the non-interacting Hamiltonian Hq, 
namely 

\^(t ))i = |n), H \n) = E n \n) . (A.10) 

To highlight the main ideas of the DRG, we will focus on the evolution of the (persistence) 
amplitude, 

C n (t) = (r#(t))i , (A.ll) 

though the time evolution of the off-diagonal overlap coefficients C m (t) = (m\ip(t))i;m 7^ n 
can be studied along the same lines. In this section we restore the diagonal matrix ele- 
ments of the perturbation allowing for (n\Hi\n) 7^ to highlight the consistency of the DRG 
resummation. 

Inserting the identity \ m)(m\ = 1 appropriately, we find 



C n (t) = (n\U(t,t )ij)(t ))i = C n (io) 

rt rt' 



A 2 

1 - i A (t - t ) (n\Hi\n) - — (t - t ) 2 (n\Hj\n) 2 



-A 2 



[ [ Yl \( n \ H i\m)\ 2 e*^-^)^"*") dt> dt" + 0(X 3 
J to J t m ^ n 



(A.12) 



where we have written C n (to) = (n\ip(to))i despite our choice of initial state for which C n (to) = 
1, to emphasize that the initial amplitude factors out. 
Introducing the spectral density 

= £ |(n|# 7 |m)| 2 S(u - E m ) (A.13) 

as in the main text, the straightforward integration over the time variables t',t" in eq.(A.12) 
can be performed, yielding 
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C n (t) = C n (to)l l-iX(t-t ) (n|flj|n> - — (t - t ) 2 (n|#,|n) 2 



-iX (t - t ) / dw 



-A 2 / du, Pn{u) 



{E n ~ U) 



sin[(.E n -u)(t- t Q )] 
(E n - u)(t - to) 



(A.14) 
(A.15) 



1 - cos[{E n - oj)(t - to)} 



+ 0(X 3 



Using eqs. (2.41,2.42), we find the asymptotic long time limit, 

r A 2 

C n (t) = C n {to)l 1 — £ A (i — t ) (n\Hj\n) - — (t - t ) 2 (n\Hj\n) 2 



n M -"m 



-iit - to) X 2 £ " A 2 vr it - to) p n (En) 

-A 2 V' 



(n\Hj\m)\^ | ( 1 3 



(E n — E m ) 2 \t — to 



(A.16) 



where J2 m ' refers to the sum over all the states with E m ^ E n . We write the expression 
above as, 



C n {t) — C n {to) 



1 + A S«(t) + A 2 S n 2 \t) - A 2 4 2) + O (j±- t A 3 



(A.17) 



where X^S^\t) are secular terms (linear and quadratic in time respectively) while the term 
—X 2 Z n 2 ^ in eq.(A.17) is the time independent term in eq.(A.16). 

To use the dynamical renormalization group, we renormalize the amplitude as 

Cn(to) = C n (T) Rn(r), R n {r) = 1 + A r« (r) + A 2 r® (r) + 0(A 3 ) . (A.18) 

Hence, 

Cn{t) = C n {r) {l + x\r n l \r)+S n l \t) \ + A 2 '\r n 2 \r) + S n 2 \t) + r«(r) 5«(t) 



-A 2 ^ 2) + 0(A a 



(A.19) 



The counterterms (r) are chosen so as to cancel the secular terms at a time scale t = r, 
leading to 



r«(r) = -SW(r) , r«(r) = -^(r) + [rW(r)] s 



(A.20) 
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The independence of the amplitude C n (t) on the arbitrary time scale r, namely dC n {t)/dr = 
leads to the dynamical renormalization group equation 



C n (r) 



C n {r) 



+ 



\rW(r) + A 2 (r< 2 >(r) + r^(r) + 0(A 3 ) 



(A.21) 



where the dot stands for derivative with respect to r. Keeping terms up to second order in A 
yields, 



C n (r) 



^{ArW(r) + A 2 [r^(r)-\ (r«(r) 



+ C(A 3 



Notice that the t-dependent pieces cancelled out identically. Such cancellation is necessary 
for the consistency of the DRG to second order in A . We find using the counterterms given 
by eq.(A.20), 

C n (r) = CM e ^W* a {«£VH . (A.22) 

Using now the explicit form of the secular terms which are read off from eq.(A.16), we finally 
find the DRG solution to second order to be given by 



C n (r) = C n (0) e" 4 



AE n 



e 2 



with the energy shift and width up to 0(A 2 ) given by 



AE n = A (n\Hi\n) + A 2 ^ 



(n\Hj\m)\' 

E — E 
J-'n J - / m 



, T n = 2tt p n (E n ) . 



(A.23) 



(A.24) 



Notice that the terms in r 2 cancel identically in the exponent of eq.(A.22) and (A.23). Again, 
such cancellation is necessary for the consistency of the DRG [18]. 

Finally, inserting eq.(A.23) into eq.(A.19) and choosing the arbitrary renormalization 
scale r to coincide with the time t, we find the asymptotic long time behavior of the dynamical 
renormalization group resummed amplitude as 



C n {t) 



t^oo „ -iAE„t Fnt 



e 2 . 



(A.25) 



The energy shift AE n is clearly the same as obtained in familiar perturbation theory, the 
decay rate T n coincides with the transition probability per unit time obtained from Fermi's 
Golden rule, and the wave function renormalization 



Z n = l- A 2 



\(n\Hi\m)\ 2 d 
(E n — E m ) 2 dE r 



(E n + AE n ) 



(A.26) 



is the same as obtained in usual perturbation theory and determines the overlap between 
the unperturbed state and the exact eigenstate of the Hamiltonian. Thus, we see that the 
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dynamical renormalization group leads to an alternative formulation of the construction of 
the exact eigenstates which allows to extract the energy shifts, widths and weights and that 
coincides with the familiar setting of quantum mechanics. Its results also coincide with those 
obtained in the Markovian approximation to the WW equations considered in the main text. 

This exercise has an important byproduct. By considering the perturbative expansion 
up to second order in the interaction Hamiltonian, we found secular terms that contain the 
square of the first order secular terms. These must cancel consistently and systematically 
in the final DRG equation since they are accounted for in the resummation furnished by 
the DRG. Indeed we found such cancellation in eqs.(A.22,A.23). This is the equivalent of 
the renormalizability of the perturbative expansion. Thus, this example not only provides a 
pedagogical framework to explore and confirm the dynamical renormalization group, but also 
manifestly shows the renormalizability in the sense that higher order terms that appear in 
the perturbative expansion, which are associated with the expansion of first order terms in 
the solution, cancel systematically in the dynamical renormalization group equation. 

B Perturbative cancelation of vacuum diagrams 

We saw in the Minkowski case that the C n 's could be renormalized in a way so as to absorb 
the effect of disconnected vacuum diagrams. In that situation this was simple to do since 
the wavefunction renormalizations were constant in time. However, this does not hold in a 
cosmological setting. Thus, the question arises as to how to deal with the effects of vacuum 
diagrams. In this subsection, we give a prescription for handling this issue. 
First, add the following counterterm to the Hamiltonian 



where Ao(r?) is given by eq.(4.33). Next, follow the Minkowski procedure by redefining the 
amplitudes as 



H c t 







iHct(ri) + i (r?) C K (rj) - i{n\Hi( V )\A) C A {v) 



(B.3) 




(B.4) 



Since Ca ~ 1 and C K oc A and because H c t, zq cx A 2 it follows that 



i{ K \Hi(i])\A} Ca(v) oc A 



while 



iH ct (ri) + z (rj) C K {rj) oc A 3 . 
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Therefore the latter terms in eq.(B.3) are perturbatively small, so that to leading order 
eq.(B.3) becomes: 

C K (rj) = -i(K\Hi(r])\A) C A (v) ■ (B.5) 

However, in (B.4) both terms are of the same order and thus we must keep both of them for 
consistency. Following the steps presented in the previous sections (after neglecting the first 
term in the right hand side of (B.3)), implementing the Markovian approximation, we finally 
obtain the equation for the single particle amplitude: 

C l0 ( V ) + [iH ct ( V ) + zo{v) + Wo,o] C lp M + Wb^fa) = 0, (B.6) 

where Wo,o is the vacuum contribution and Wo,i is the single particle contribution eq.(4.43). 
It is clear that the vacuum contribution is now canceled since 



H ct {rf) = -Im[W ,o], iofo) = -Re[W ,o] • (B.7) 

The solution of eq.(B.6) with eq.(B.7) is given by eq.(4.47). It should be clear that the same 
procedure cancels the vacuum contribution in the equation of motion for the renormalized 
amplitudes for all excited states. 
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